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ABSTRACT 

Large-scale modes in the temperature anisotropy power spectrum C/ measured by 
the Wilkinson Microwave Anisotropy Probe (WMAP), seem to have lower amplitudes 
(C2, C3 and C4) than that expected in the so called concordance ACDM model. In 
particular, the quadrupole C2 is reported to have a smaller value than allowed by 
cosmic variance. This has been interpreted as a possible indication of new physics. 
In this paper we re-analyse the WMAP data using the 2-point angular correlation 
and its higher-order moments. This method, which requires a full covariance analysis, 
is more direct and provides better sampling of the largest modes than the standard 
harmonic decomposition. We show that the WMAP data is in good agreement (~ 30% 
probability) with a ACDM model when the WMAP data is considered as a particular 
realization drawn from a set of realistic ACDM simulations with the corresponding 
covariance. This is also true for the higher-order moments, shown here up to 6th order, 
which are consistent with the Gaussian hypothesis. The sky mask plays a major role 
in assessing the significance of these agreements. We recover the best fit model for the 
low-order multipoles based on the 2-point correlation with different assumptions for 
the covariance. Assuming that the observations are a fair sample of the true model, 
we find C2 = 123 ± 233, C3 = 217 ± 241 and C4 = 212 ± 162 (in /iif^). The errors 
increase by about a factor of 5 if we assume the ACDM model. If we exclude the 
Galactic plane \h\ < 30 from our analysis, we recover very similar values within the 
errors (ie C2 = 172, C3 = 89, C4 — 129). This indicates that the Galactic plane is not 
responsible for the lack of large-scale power in the WMAP data. 

Key words: cosmology - cosmic microwave background 



1 INTRODUCTION 

Measurements of the cosmic microwave background (CMB) 
anisotropies by WMAP (Bennett et al. 2003; Spergel et al. 
2003; Peiris et al. 2003) are in good agreement with a 'con- 
cordance' cosmology based on the ACDM model. However, 
the WMAP results also confirm the low amplitude of the 
CMB quadrupole first measured by COBE (eg Hinshaw et 
al. 1996a). As Bennett et al. (2003) comment, the ampli- 
tudes of the quadrupole and the octopole are low compared 
with the predictions of ACDM models. The WMAP team 
also present a convincing case that the low CMB multi- 
poles are not significantly affected by foreground Calactic 
emission (see also Tegmark, de Oliveira-Costa and Hamil- 
ton 2003). The discrepancy between the observations and 
the ACDM model is particularly evident in the tempera- 
ture angular correlation function, which shows an almost 
complete lack of signal on angular scales ^ 60 degrees. Ac- 



cording to Spergel et al. (2003), the probability of finding 
such a result in a spatially-fiat ACDM cosmology is about 
1.5 X 10"^. This is small enough to require an explan ation, 
e.g. in the form of new ph ysics (IContaldi et al. 20ol or m 
terms of spatial curvature llEfstathioi^003all . 



In this paper we consider the WMAP data from the 
point of view of the two-point angular correlation function 
as opposed to the usual harmonic decomposition. In some 
aspects our approach is similar to the study presented by 
Hinshaw et al. (1996b) for the COBE data. The paper is 
organized as follows. In §2 we introduce our methodology 
and test it with realistic simulations. In §3 we apply our 
method to the real WMAP data and present a comparison 
with the simulations. We also describe the application of 
higher-order statistics to test the Gaussianity of WMAP. 
The conclusions are described in §4. 
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2 2-PT CORRELATION W2{e) IN 
ACDM MODELS 

On the largest scales, the 2-point correlation W2{6) has some 
important advantages and disadvantages over spherical har- 
monics. The overall shape of the 2-point function is very sen- 
sitive to small differences in power at low multipoles. This 
is clearly an advantage when measuring the quadrupole. On 
the other hand, different bins in 6 are highly correlated, 
which means that we need to use the full covariance matrix 
to assess the significance of any departures between measure- 
ments and models. Another important advantage of using 
the 2-point function is that large-scale modes can be easily 
sampled from any region of the sky, even when a mask is 
required, provided the region is large enough. In compari- 
son, harmonic decomposition breaks the angular symmetry 
in the sky into orthogonal bases which have an arbitrary 
phase orientation. This is not a problem in the case of full 
sky coverage, but the results will be coordinate-dependent if 
a mask is used. Moreover, the harmonic coefficients aim (see 
Eq|^ over masked data are convolved with the harmonic im- 
age of the mask. The resulting coupling matrix needs to be 
inverted and the resulting estimator for the recovered coeffi- 
cients is typically biased because of the non-linear transfor- 
mations involved. These potential complications are avoided 
in configuration space. 



2.1 Definition and estimators 

The 2-point angular correlation function is defined as the 
expectation value or mean cross-correlation of density fluc- 
tuations at two positions qj^ and q2 in the sky: 

W2ie) = (Ar(qi)Ar(q2)), (i) 

where 6 = |q2 — qi|, assuming that the distribution is sta- 
tistically isotropic. 

To estimate 102 (9) from the pixel maps we use a simple 
but effective estimator: 



Table 1. CMBFast Cosmological Parameters 



W2{e) = 



(2) 



where Ai and Aj are the temperature differences in the map 
and the sum extends to all pairs separated by ^ ± A9. The 
differences are normalised so that (Ai) = 0. The weights Wi 
can be used to minimise the variance when the pixel noise 
is not uniform, but this introduces larger cosmic variance. 
Here we follow the WMAP team and use uniform weights 
(i.e. Wi = 1). We use bins /S.9 whose size is proportional to 
the square root of the angle 6. 

This estimator is unbiased and is equivalent to the min- 
imum variance estimator of Landy & Szalay (1993). For 
higher-order moments see §3.6. 



2.2 The ACDM Model and simulations 

Temperature fluctuations can be expanded in terms of spher- 
ical harmonics according to. 
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72 km/s/Mpc 




.046 




.224 


Ha 


.730 




0.0 


T 


.17 


Vs 


1.0 



which deflnes the multipole power spectrum: Ci =< 
\aim\^ >• Gaussian simulations can then be made by taking 
aim to be independent Gaussian flelds with random phases, 
zero mean and variance C;. 

In order to test the signiflcance of the correlations mea- 
sured in the WMAP data, we make simulated sky maps by 
using the best-fit cosmological parameters estimated from 
WMAP (Spergel et al. 2003). The theoretical power spectra 
are generated using the publicly available CMBFAST pack- 
age (Seljak & Zaldarriaga 1996). The CiS are normalised to 
the COBE spectra and the input parameters are listed in Ta- 
ble Q Along with cosmological parameters, CMBFAST ac- 
cepts a reionization optical depth, r, and the scalar spectral 
index, 7)^ . From now on we refer to this particular set of para- 
menters as the ACDM model. Figure Q shows the resulting 
power spectrum for the ACDM model. In the ACDM model 
the first three non-zero multipoles are (in f^K^): 

C2 = 1130, C3 = 537, C4 = 304, (4) 

while Bennett et al. (2003) find from WMAP data: 



C2 = 129 ± 799 
C3 = 320 ± 320 
C4 = 316 ± 104. 



(5) 



AT(q) 



EE 

t = m=-i 



(3) 



This power spectrum is then used as input to generate 
all-sky CMB maps using the HEALPix^ package (Gorski, 
Hivon & Wandelt 1999). We simulate 100 V-band maps 
(FWHM = 21' at 61 GHz), each with a different random 
phase, using the best fit WMAP power spectra^ These maps 
are all generated with the RING pixelization scheme and 
nside=512. Noise is added to the simulations following the 
pixel-noise pattern in the WMAP data. We also impose the 
most conservative WMAP kpO foreground mask to the sim- 
ulations. 

For some crucial aspects of the analysis we have used an 
additional set of 2000 V-band ACDM simulations to check 
how the results depend on the number of simulations used. 
In general we find that with 20 simulations one finds similar 
qualitative conclusions to those with 2000 simulations, with 
a typical uncertainty of ~ 50% in the intervals of confidence 
(see below §3.4). 



^ http:/ /www. eso.org/science/healpix/ 

^ We have also generated 20 W-band simulations to compare with 
W-band WMAP, but find no significance differences with V-band 
and decided to use only the latter because of the smaller pixel 
noise. 




Figure 1. Power spectrum C; for the ACDM model (continuous 
lino). Open squares show the first multipoles by Bennett et al. 
(2003). Closed triangles with errors correspond to Eg. llMI . Unlike 
the conventional way of plotting this curve, each amplitude here 
shows the contribution of each multipole to the sky anisotropics 
in Eq. |7|. Dashed line shows the low-Q model in Eq. 

2.3 The low-Q ACDM model 

As a useful comparison model, we design a fiducial variation 
of the ACDM model with low values of the quadrupole and 
octopole. This new model is identical to the ACDM model 
for I > 3 but has zero quadrupole C2 ~ and half the 
octopole C3 — > C3/2. As we will show later, this model 
is in better agreement with the WMAP data than the 
ACDM model an d follows the curv ed universe with a trun- 
cated spectrum iEfstathiou 2003al) . Thus, we modify the 
ACDM power spectrum by setting: 



^low- 
^2 


Q 


= 
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1 ^ACDM 

- 2^' 




Q 


= C{'^''", I > 3 



We refer to this model as low-Q ACDM . These new Cj's 
are then used to generate 20 HEALPix V-band maps, each 
with a different random phases. 

2.4 Predicted W2{e) from C; 

The 2-pt function, W2{6), is the Fourier transform of the an- 
gular power spectrum of temperature fluctuations. In terms 
of the Ci in Eq. 0, we have (see e.g. Bond & Efstathiou 
1987): 

W2{e) = ^^ {2l + l)CiPi{cose), (7) 

where Pi are the Legendre polynomials. The sum starts at 
I = 2 because the monopole and the dipole have been sub- 
tracted from the WMAP data (CMB rest frame). 

Figure shows the multipole coefficients for the 
ACDM and low-Q models. Figure|21shows the corresponding 



Figure 2. Theoretical prediction of the 2-point angular corre- 
lation W2 in the ACDM model with all multipoles (continuous 
line), without the quadrupole (long dashed line) and without 
quadrupole and octopole (short dashed line). 



analytical summation for all multipoles I ^ 2 (continuous 
line), / 3 (long-dashed line) and I ^ 4 (short-dashed line). 
The figure very clearly illustrates the sensitivity of the angu- 
lar 2-point function to the low multipoles: e.g. varying just 
the quadrupole completely changes the shape of the curve 
at large scales. 



2.5 •W2{9) from simulations 

Figure |^ compares the theoretical values of W2 (6) in the 
ACDM and low-Q models with that estimated from the sim- 
ulated maps. The continuous line is the mean of the different 
realizations of each model and the dashed line is the theo- 
retical estimate corresponding to the set of Cjs. The shaded 
area shows the 1-sigma dispersion of the simulations in each 
bin. As can be seen from the figure, the ACDM model has 
the larger errors. This is expected since the ACDM model 
has stronger lower-order multipoles which in turn lead to 
larger variations from realization to realization. 

The right-hand panels in Figure |3 show the estimation 
of W2 for pixels within Galactic latitudes |&| > 30 degrees. 
This reduces the number of pixels by almost a factor of two 
which results in larger errors. Note how the recovered values 
of W2 are not strongly biased within the error bars. Thus, 
even after excluding a broad Galactic region one can still 
use W2 to separate a low and high quadrupole in the data. 
We will further elaborate on this point later in the paper. 



2.6 Error bars 

The absolute error (1-sigma dispersion) as a function of scale 
from different realizations of each model is shown in Figure 
13] As can be seen from the figure, errors for the ACDM simu- 
lations (long-dashed lines) are much larger than errors in the 
low-Q model (short-dashed line). Note also how the low-Q 
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Figure 3. The 2-pt function W2{d) from simulations (continu- 
ous line) with 1-sigma confidence region (shaded region) com- 
pared to W2{d) from the input C; spectrum (dashed line). Top 
left: ACDM model with kpO-mask. Top right: ACDM model with 
kpO-mask and |fe| > 30 Galactic cut. Bottom left: low-Q model 
with kpO-mask. Bottom right: low-Q model with kpO-mask and 
|6| > 30 Galactic cut. 



model estimation is more noisy since it is based on 20 re- 
alizations, while there are 100 and 1000 realizations in the 
ACDM model. 

The shaded region represents the dispersion (2-sigma) 
in the jackknife errors (see below) for the different low-Q 
realizations. 



2.7 Covariance Matrix 

As mentioned before there is strong covariance between dif- 
ferent bins in W2(^). This is mostly due to the fact that the 
lowest I multipoles, and therefore larger scales, have more 
power than higher I multipoles (e.g. see Fig.0. As W2{Q) is 
dominated by large scale modes, the correlations at different 
scales are strongly correlated. This is obviously worse in the 
ACDM model than in the low-Q ACDM , which has a lower 
quadrupole and octopole. 

It is therefore essential that we estimate the covariance 
between different bins in ■W2{0)- This can be calculated from 
the simulations by using the following definition of the co- 
variance matrix: 



Cij = {Aw2{e,) Aw2ie,)) 

N 

= ^$]At«2^(e,)A^2^( 

L=l 

where 

Aw^ie,) = w^{e,)-^2(ei). 



(8) 



(9) 




Figure 4. Short, long-dashed and dotted lines show the absolute 
errors (or dispersion) in W2 from the low-Q simulations, the first 
40 and the first 1000 ACDM simulations respectively. The con- 
tinuous lines corresponds to estimated jackknife errors with N=8 
and N=32 subsamples within the WMAP data, while the shaded 
region shows the dispersion in jackknife errors from each of the 
low-Q realizations. 



the A'^ realizations. The case i = j gives the error variance: 
(Tj = di. 

It is also useful to define the normalised covariance ma- 
trix (or correlation matrix), as: 



Cij — 



(10) 



Here W2 {Oi) is the 2-point function measured in the L-th 
realization [L — 1 . . . N) and W2{9i) is the mean value for 



•\/ CiiCjj 

Figure |K| shows the normalised covariance in the 
ACDM (left) and low-Q model (middle). The covariance 
shown for the ACDM model corresponds to 40 realizations, 
but very similar pictures are obtained for a 1000 realizations. 
The covariance matrix is quite different in these two cases 
as the ACDM model shows correlations on larger scales be- 
cause of the higher quadrupole and octopole. Note in partic- 
ular the strong anticorrelation between angular bins in the 
range 50-100 degrees with the bins on larger and smaller 
scales. This is clearly an indication of very large variations 
in the low order multipoles within different realizations (eg 
see Fig|5J. Note also how the errors in the first bin (zero- 
lag in W2) is dominated by instrumental noise so that this 
bin is not correlated with the rest of the bins (grey color 
corresponds to zero). 

2.8 Jackknife covariance 

We can also estimate the errors in W2{0) from a single sim- 
ulation, or from the real sky with a variation of the jack- 
knife error scheme proposed by Scranton et al. (2001). This 
has the potential advantage of producing an error estimate 
which is model independent. In the estimation, the sample 
is first divided into M separate regions on the sky, each 
of equal area. The analysis is then performed M times, on 
each occasion removing a different region. These are called 
the (jackknife) subsamples, which we label k — 1 . . . M. The 
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Figure 5. Normalised covariance (or correlation) matrices between different angular bins in W2{d) . Grey scale goes from -1 (black) to 
+1 (white). Left: ACDM model. Middle: low-Q model. Right: jackknife samples in low-Q model. 



estimated statistical covariance for W2 at scales 6i and 6j is 
then given by: 



M ■ 



M 

1 



M 



■ W2{ 



(11) 

(12) 



where W2{0i) is the measure in the fc-th subsample (fe = 
1 . . . M) and W2{6i) is the mean value for the M subsam- 
ples. The case i — j gives the error variance. Note how, 
if we increase the number of regions M, the jackknife sub- 
samples are larger and each term in the sum is smaller. We 
typically take M = 16 corresponding to dividing the sphere 
into 8 octants and each octant in two (by longitude). We 
have checked that the resulting covariance gives a stable an- 
swer for different choices of shapes and M from M = 8 to 
M = 32. 

The mean normalised jackknife covariance shown in the 
right panel of Fig. |S] is in very good agreement with the 
sample covariance shown in the middle panel. The jackknife 
covariance appears smoother because it is the mean of 20x16 
jackknife samples, while the sample covariance is the mean 
of only 20 samples. Other than this, the jackknife matrix 
correctly captures all the relevant correlation information 
present in the ensemble. The corresponding diagonal errors 
are shown in Fig. |1] 

In general, one would not expect such a good agree- 
ment, because the jackknife method can not account for 
variations on scales larger than the jackknife samples. In 
our case, however, by construction there is little power on 
the largest scales (quadrupole and octopole) in the low-Q 
model. This explains the good performance of the jackknife 
method. 



2.9 X test 

In order to quantify how well a model fits the data we use a 
test. Since the bins of the W2{d) are correlated, however, 
the simple x^-test is not valid but must be modified by using 
the covariance matrix. The value of the x^ is given by 



(13) 



where Ai = (Si) — ■W2' (6i) is the difference between the 
"observation" O and the model M. In terms of the reduced 
covariance matrix, this is 



2 

X = 



E (I) m 



(14) 



2.10 Singular Value Decomposition 

In using real world data, one must worry about degeneracies 
in the covariance matrix due to an over-determined system 
and noise. If the rows (or columns) of dj are not indepen- 
dent, the determinant of the matrix is zero and there is no 
inverse. With real data, the determinant of dj will not typi- 
cally be zero but the matrix can have very small eigenvalues 
which then lead to artificially large eigenvalues of the inverse 
matrix. 

In order to eliminate the degeneracies in the covariance 
matrix, we perform a Singular Value Decomposition (SVD) 
of the matrix, 



a, = {u,k)^WkiVij, 



(15) 



where Wki is a diagonal matrix with the singular values on 
the diagonal. By doing the decomposition, we can choose the 
number of modes we wish to include in our x^ by effectively 
setting the corresponding inverses of the small singular val- 
ues to zero. 

In practice we find that the SVD is only required when 
we use a small number of simulations to estimate the covari- 
ance matrix. With more than 100 simulations we find that 
the covariance matrix inversion is stable and one then gets 
more accurate x^ values. 



2.11 Recovering Ci from ■W2{0) 

Using the SVD of the covariance matrix we can recover C;s 
from W2{6) by requiring a minimum x^ (see Szapudi et al. 
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Figure 6. Recovered values of (C2,C3) (upper panel) and 
(C2,C4) (lower panel) from the ACDM model. The black box 
indicates the values used to generate W2{6) and the gray box the 
recovered value. Confidence contours are plotted at 10, 68 and 
98% of the minimum value. 
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Figure 7. Recovered values of (C2,C3) (upper panel) and 
{C2,C4) (lower panel) from the low-Q model. The black box in- 
dicates the values used to generate W2{0) and the gray box the 
recovered value. Confidence contours are plotted at 10, 68 and 
98% of the minimum value. 



2001 for a different approach). To test the method, we have 
recovered the low C;S from the ACDM and low-Q modes. 
These are shown in Figures Fig. |n] and Fig. Q respectively. 
The confidence contours are plotted relative to the mini- 
mum value. As it is clear from the figures, the method 
accurately reproduces the correct values of the C;s. We have 
also checked that this method is robust and insensitive to 
the number of modes used in the calculation. 



3 RESULTS FROM WMAP 

3.1 Estimation of ^12(6') in WMAP 

Figure |H| shows different estimations of W2{0) from the 
WMAP data. These different estimates gives an idea of 
the systematics involved, which are related mostly to pos- 
sible Galactic and foreground contamination in the maps 
(see Bennet et al. 2003, Tegmark et al. 2003 for more de- 
tails). Also note how the jackknife errors (from the disper- 
sion in the V-band estimates from different regions in the 
sky) roughly include the variations in these different esti- 
mations. Thus we will use V-band estimates with jackknife 
errors as representative of WMAP sampling variance and 
systematic uncertanties. 

Note how our estimates of W2(0) agree well with that 
presented by the WMAP collaboration (e.g. compare to 
Fig. 13 in Bennett et al. 2003). 
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Figure 8. Error bars represent the 2-point angular correlation 
and 1-sigma jackknife dispersion W2 in the WMAP V-band. The 
continuous line, passing through the error-bars, shows U12 in the 
WMAP W-band. These two cases have the kpO mask, while the 
dotted line excludes all the Galactic plane with b > 30 deg. The 
long and short dashed lines correspond to the foreground clean 
map of Tegmark etal (2003) with and without a Galactic cut 
\h\ > 20. 



3.2 Comparison of W2{6) to simulations 

Figure |n| compares the WMAP W2{9) estimation with the 1- 
sigma dispersion in the corresponding ACDM simulations. 
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Figure 9. The 2-point angular correlation W2 in the WMAP 
V-band data (points with jackknife error bars) compared to the 
ACDM simulations (shaded region corresponds to 1-sigma error 
dispersion). In both cases we use the kpO Galactic mask. 

At first sight one might conclude that a ACDM model is 
ruled out to a high significance. It should be remembered, 
however, that there is a strong covariance in the W2{0) es- 
timation (see Fig|3 . This point is illustrated in Figure 1101 
which shows 10 of the ACDM realizations that happens to 
have a low quadrupole and are compared with one that hap- 
pen to have a high quadrupole (all chosen from the first 
100 ACDM simulations). Note how systematics effects in 
WMAP (which are well represented by the jackknife er- 
rors) are comparable to the differences between these 10 
ACDM simulations and the mean WMAP estimation. 

These very obvious differences in the lowest multipoles 
in some of the realizations are also apparent in the all-sky 
simulated maps. Figure 1111 illustrates this point by show- 
ing one of the low quadrupole ACDM realizations and a 
ACDM realization with the average value of the quadrupole. 

In Figure [T^ we show the comparison of W2{6) with the 
low-Q model which is apparently in much better agreement 
with WMAP. 

3.3 Covariance Matrix 

To quantify these differences more accurately we need to 
calculate the covariance matrices of the models. We can 
also estimate the covariance matrix of the data using the 
jackknife method, Eq. 1111 . Figure [T^ shows the normalised 
covariance matrix estimation, which is more similar to the 
low-Q model than to the ACDM model shown in Fig.|S| Fig. 
|1| compares the diagonal errors with N=8 and N=32 jack- 
knifes subsamples (continuous lines) with the errors in the 
low-q and ACDM models. 

3.4 tests of tlie WMAP data 

Using the method described in the previous section, we com- 
pare the WMAP W2{e) to different models in Tabled The 
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Figure 10. The 2-point angular correlation U12 in the WMAP 
V-band data (points with jackknife error bars) compared to 10 of 
the ACDM simulations (continuous lines) which happens to have 
a low quadrupole and one (dashed line) with a high quadrupole. 
These are all taken out of the 100 realisations, whose dispersion 
is shown as the shaded region in Figlol 




Figure 11. Two of the ACDM simulations, with average (top) 
and low (bottom) quadrupole. 



second row shows how considering the WMAP data as a 
realization of the ACDM model gives a very good fit (62% 
C.L.). However, considering the WMAP jackknife errors'^ 



A similar result is obtained using low-Q ACDM true errors. 
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Figure 12. The 2-point angular correlation W2 in the WMAP 
V-band data (points with jackknife error bars) compared to the 
low-Q simulations (ACDM model with zero quadrupole, C2 = 
0, and 1/2 the octopole C3/2). The continuous line shows the 
ACDM model predictions. 
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Figure 13. Normalised covariance matrix from WMAP jackknife 
samples, to be compared to Fig. l5l 



Table 2. Likelihood of WMAP data in different models. First col- 
umn shows the model we are comparing to, e.g. ACDM denotes 
that we are comparing to the mean W2{6) of the ACDM simula- 
tions. The second column indicates the model from which we take 
the covariance matrix. The fourth column indicates how many 
singular values we use in calculating the ^^nd the fifth column 
is the probability of the fit. 



Model 


Errors 




Modes 


P 


ACDM 


WMAP 


200 


7 





ACDM 


ACDM 


3.5 


5 


0.62 


low-Q 


WMAP 


9.6 


7 


0.21 


low-Q 


low-Q (true) 


4.7 


5 


0.45 


low-Q 


low-Q (jackknife) 


4.1 


6 


0.66 



Table 3. Likelihood of WMAP data with a galaxy cut in different 
models 



Model 


Errors 




Modes 


P 


ACDM 30 


WMAP 30 


191 


7 





ACDM 30 


ACDM 30 


5.7 


6 


0.46 


low-Q 30 


WMAP 30 


14 


7 


0.06 


low-Q 30 


low-Q 30 


5.2 


7 


0.64 



(i.e. first row in Table 2) leads to a vanishing probability for 
the ACDM model, as was already hinted at in Figl^ 

On the other hand, the WMAP data provides a res- 
onable fit to the low-Q models, even if we use the WMAP 
errors, as can be seen from the third, fourth and fifth row 
of Table 13 The WMAP data clearly has a low quadrupole 
compared to ACDM as was already suggested by Figure 13. 

Considering only the WMAP data at \b\ > 30 deg, and 
comparing the resulting W2{9) to the models, give similar 
results as can be seen in Table|3 This is in qualitative agree- 
ment with that found by COBE (Kogut et al. 1996). 

The above results are obtained with a relatively small 
number of simulations (40 for ACDM and 20 for low-Q 
model) and also a small number of jackknife subsamples 
(16). Thus, although these estimates probably give the right 
order of magnitude for the calculation we can not expect 
them to be very accurate. We next quantify what are the 
errors involved by repeating the ACDM analysis with an in- 
creasing number of simulations, from A'^ = 100 to A'^ = 2000, 
to calculate the covariance matrix. We can also check with 
these additional simulations if the distribution gives an 
accurate representation of the probabilities. This is illus- 
trated in Fig ll4l which shows the histogram of values ob- 
tained for each ACDM realization using a direct inversion of 
the covariance matrix (no SVD), estimated from A'" = 1000 
ACDM realizations. The histogram matches well the proba- 
bilities given by the x^ distribution. Moreover we find that, 
other than Poisson errors, the values are quite insensitive 
to the total number TV of simulations used. The probability 
of finding a ACDM simulation with a x^ value larger than 
WMAP seems to converge to 32%. This number compares 
well to the crude 62% estimate we found in Table H using 
SVD over 40 simulations. If we used 20 or 100 simulations 
with SVD we find intermediate results. We have also esti- 
mated the x^ values for A'^ = 2000 ACDM simulations using 
the WMAP jacknife covariance matrix. We find that none 
of the A*' — 2000 ACDM simulations have a x^ smaller than 
WMAP (the smaller value in the 2000 simulations is x^ = 62 
for 36 bins in w^iO) ). Thus the probability for ACDM given 
the WMAP errors is P < 0.05%, in good agreement with 
first entry in Table |21 



3.5 Recovered Ci from W2{p) 

Using the previously discussed method (see §2.11) of calcu- 
lating the with the SVD, we can recover the CiS from the 
data by requiring that x^ has a minimum. Fitting three mul- 
tipoles (C2, C3 and C4) of a ACDM model to the WMAP 
data and using the ACDM errors yields 



C2 = -7 C3 = 190 C4 = 226. 



(16) 
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[for 36 angular bins in vr^iO) ] 

Figure 14. Number of ACDM simulations with given values of 
in Eg 1131 using the covariance matrix from the same ACDM sim- 
ulations. The same test for WMAP gives: — 39 (shown by the 
arrow), which is below 32% of the simulations. 
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Figure 15. Confidence plot of (C2,C3) (upper panel) and 
(C2,C4) (lower panel) recovered from the ACDM model with 
WMAP (inner contours) and ACDM (outer contours) errors. The 
contours, in both cases, are at the 0.1, 0.68 and 0.98 confidence 
levels. The filled square shows the location of the theoretical 
ACDM values in both cases. 



With the jacknife WMAP errors we get 
Ca = 123 Cg = 217 d = 212. 



(17) 



Comparing with Eq. 0] we see that the best fit value has a 
negligible quadrupole and an octopole that is roughly half 
the value of the ACDM model. Even the d is somewhat 
lower than the ACDM value. The low-Q model is hence 
a better approximation to the W2{6) measured from the 
WMAP data. 
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The confidence contours of the Cis are plotted in Fig. 
1151 for the WMAP errors for the two sets of pairs (C2,C3) 
and (C2,C4). The figure fully confirms our expectations: 
the WMAP data is in good agreement with ACDM when 
WMAP data is considered a realization of the ACDM model 
(i.e. using ACDM errors). The WMAP data is not in agree- 
ment with ACDM however when using WMAP errors. 

From the confidence plots we can recover the 68% C. 
L. (or ~ 1-sigma) limits for the C;s by considering the pro- 
jections of the contours on the axes. When no model is as- 
sumed, i.e. we use WMAP errors, we find (in fiK^): 



C2 
C4 



123 ± 233 
217 ± 241 
212 ± 162, 



(18) 



which are shown as closed triangles in Fig0 With a 
ACDM covariance matrix the errors are much larger (in 



C2 = -7 ± 1568 
C3 = 190 ± 799 
C4 = 226 ± 339, 



(19) 



Compared to Eq|K| our errors are about a factor of two 
larger. 

3.6 Higher-order moments 

We also present the higher-order moments of the 2-point 
correlations in WMAP. This provides a consistency test for 
the Gaussian hypothesis, which is implied in the comparison 
with simulations. 

The higher-order moments of the 2-point function can 
be defined as: 



kp,{e) = (Ar''(ei)Ar«(e2)) 



(20) 



where p = q = 1 corresponds to the 2-point function in 
Eq.0, i.e. k\\{9) — W2{0) and p = gives the 1-point mo- 
ments, e.g. ko2 = fcii(O) is the variance and fco3 = fci2(0) is 
the skewness. It is then useful to define 2-point cumulants 
as connected moments, subtracting the lower-order contri- 
butions (e.g. see Gaztanaga, Fosalba & Croft 2002) : 



Wi2{0) 

wri{0) 



W2{e) ; 

ki2{e) 
ki3{e) 
ki4{e) 



fcii(e') 

-3fcii(o) kii{e) 

-6fcii(0) fci2(e) -4fci2(0) fcii(e) 



(21) 
(22) 
(23) 
(24) 



For a Gaussian field we expect Wpq{9) = 0. 

It is straightforward to estimate these higher-order mo- 
ments in the same way in which we have calculated the 2- 
point function: 



P Q 



(25) 



Here, as in calculating the 2-point function, uniform weights, 
Wi — 1, are used. We have estimated the higher order mo- 
ments as well as the covariance matrix both from WMAP 
and simulated data. 

Figure 1161 shows 1012, W13 and ^114 for the WMAP 
V-band data and the combined foreground clean map in 



10 Gaztanaga, Wagg, Multamdki, Montana, Hughes 




50 100 1.50 50 100 ^.50 50 1.00 '50 



^{degrees) (7(dc£2rccs) O(dcgrccs) 

Figure 16. Higher-order 2-pt moments in WMAP V-band (continuous line) and WMAP foreground clean map (long-dashed line) as 
compared with the dispersion in the corresponding measurements in the low-Q ACDM simulations (shaded regions). The short-dashed 
line in the left-hand panel shows Wi2 from the full combined clean map, ie without the Galactic cut at |b| < 20. 



Tegmark et al. (2003) with a Galactic cut excluding |6| < 20 
(long dashed line). Also shown are the 68% CL in the dis- 
persion from the low-Q ACDM simulations. Both WMAP 
and the combined map with |6| > 20 agree well with the 
Gaussian hypothesis Wpq = at high significance. 

For the 3rd order moment, Wi2, we also show the com- 
bined clean map without the Galactic cut (short-dashed 
line). This result is incompatible with the Gaussian low-Q 
ACDM simulations (and covariance matrix) on small angu- 
lar scales, 9 < 20 deg. This is in qualitative agreement to 
what Chiang et al. (2003) have recently found using phase 
correlators over the very same map. Nevertheless, we note 
how this apparent discrepancy goes away when we exclude 
the |6| < 20 Galactic region. Thus, this non-Gaussian feature 
seems to be related to residual Galaxy contamination, which 
is visually apparent in the image of the combined clean map 
(e.g. see Fig. 1 in Tegmark et al. 2003). 

Our best constraints on non- Gaussian features come 
from the 4th-order statistic, W13, because in this case the 
pixel noise at two different points would tend to cancel out, 
as in the 2nd-order moment W2. This is not the case for the 
3rd and 5th-order moments, W12 or wn, where pixel noise 
is added in quadrature. One needs to consider multi-point 
correlations to reduce the noise contribution and get better 
constraints on non-Gaussian models, which is beyond the 
scope of this paper (see Hinshaw et al. 1995, Komatsu et 
al. 2003). The point here is that the higher-order 2-point 
anisotropics in WMAP are perfectly compatible with the 
Gaussian hypothesis, given the pixel noise. We find similar 
results up to 8th-order moments. 



4 CONCLUSIONS 

In this paper we have studied the lowest-order multipoles in 
the CMB from recent WMAP data. This work is motivated 
by the fact that the low order multipoles, and especially 
the quadrupole as measured in harmonic space, show less 
power than that expected in the standard ACDM model 
(e.g. see Bennett et al. 2003; Spergel et al. 2003; Peiris et al. 
2003). Instead of looking at the power spectrum, we study 
the WMAP data by using the two-point angular correlation 
function, W2{9) , as a tool to estimate the multipoles. The 



main advantage is that W2 (6) is very sensitive to the low or- 
der multipoles The main disadvantage is that different bins 
of W2{0) are highly correlated so a covariance analysis is 
required. 

We have simulated a set of V-band maps with differ- 
ent random phases using the best-fit WMAP power spec- 
tra and pixel errors, and calculate the 2-point function for 
each model. Using this set of simulations we calculate the 
mean value and variance at each bin as well as the covari- 
ance matrix. The errors and the covariance matrix of the 
WMAP data is calculated by using the jacknife method. A 
ACDM model with a zero quadrupole and a low octopole is 
also simulated and used to test the reliability of the jacknife 
errors. 

Using a test with a covariance matrix we calculate 
the significance of the WMAP data in different models. It is 
found that indeed, the WMAP data with the WMAP errors 
does not fit the ACDM model and is strongly disfavoured. 
Instead, the WMAP data fits well the low-Q model with 
a vanishing quadrupole. However, taking a different point 
of view and considering the WMAP data as a sample re- 
alization of the simulated ACDM maps, we find that the 
WMAP results are fully consistent with the ACDM model. 
From this viewpoint, our sky just happens to have a low 
quadrupole and is just as likely to have a large quadrupole 
as demonstrated by the different realizations in Figure 1111 

Our results are in apparent contradiction with the har- 
monic space analysis presented by Spergel et al. (2003), 
which gives WMAP a probability as low as 0.15% of being 
a realization of the ACDM model. In contrast we estimate 
this probability to be 32% comparing W2{6) in WMAP to 
W2{9) in 1000 ACDM realizations (see Fig^. Our proba- 
bility is closer, but still larger than the 5% value obtained 
from harmonic analysis by Tegmark etal (2003)* Why are 
all these estimates so different? To answer this question we 
conduct a further test by repeating our calculation using a 
different galactic mask. We use a sharp cut in Galactic lat- 
itude |ti| of about ±20 degrees, excluding the same number 

* While this paper was being referee, we became aware of an- 
other calculation by Efstathiou (2003b) who find intermediate 
probabilities of 1 — 2%, also using direct harmonic analysis. 
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of pixels as the kpO mask. We then recalculate the covari- 
ance matrix and its inversion with this new Galactic cut 
over the same 1000 ACDM realizations. We repeat the 
analysis (in Fig |14t for WMAP and the 1000 simulations 
with this new covariance matrix. We find that only 1% of 
the simulations with this Galactic cut have a larger than 
the one measured for WMAP. Thus the final probability de- 
pends strongly on the shape of the mask, and therefore on 
how the mask is accounted for in the comparison between 
the models and observations. This is not surprising because 
the large scale modes that determine the shape of W2 {0) (eg 
FigEl extend over angular regions that are comparable to 
the mask. To recover the low order multipoles from masked 
data using harmonic decomposition, one needs to invert the 
convolution with the harmonic image of the mask. This is a 
non-linear transformation that is likely to introduce biases 
in the recovered low multipoles. We speculate that the origin 
of the above discrepancies with the harmonic analysis could 
depend on how the mask is taken into account on the differ- 
ent estimators. Because the ■W2{0) estimator is independant 
of the mask we believe that our ~ 30% probability value is 
reliable. 

One can also recover the values of the low multipoles 
from the W2{6) WMAP estimations. The validity of the 
method is demonstrated by recovering known Cts from a set 
of simulated data. Again, we find that the best fit model has 
a low quadrupole as well as a small octopole. Confidence lev- 
els in the C; plane point to a similar conclusion: the WMAP 
data fits the ACDM model, but only with ACDM priors. 
We can also reach a similar conclusion by asking how many 
of the first 100 ACDM simulations show a 2-point function 
that "looks" closer to the WMAP data than to ACDM W2 
(e.g. see Figur JTITll . There is a clear distinction of such cases, 
which is even apparent by a direct inspection of the pixel- 
map images (see Figure ITTl . We find at least 15 similar ex- 
amples of low quadrupole amongst the first 100 simulations, 
confirming our more quantitative analysis above. 

We have also considered the presence of possible non- 
Gaussianities in the WMAP data by calculating higher- 
order moments of the 2-point correlation. We find that both 
the WMAP and combined foreground clean map of Tegmark 
et al. (2003) at |6| > 20 are in excellent agreement with 
the Gaussian hypothesis. It is also shown how the full-sky 
foreground clean map (i.e. including \b\ < 20 ) shows a non- 
Gaussian feature at angular scales 6 < 20 degrees, indicating 
a possible residual galaxy contamination in the map. 

In summary, the low value of the WMAP quadrupole 
is confirmed here by looking at the two-point angular cor- 
relation function, W2 (9) , which is shown to be a well-suited 
method for the study of low multipoles in the CMB data. 
This does not necessarily indicate a need for new physics 
or non-standard cosmology, since we also show the WMAP 
W2{6) fits well a standard ACDM cosmology if one considers 
our universe as a realization of the set of possible ACDM out- 
comes. We stress nevertheless, that by considering the er- 
rors from the WMAP data alone, i.e. with the jacknife er- 
rors, 1112(6) in WMAP is strongly inconsistent with the mean 
ACDM model and, under such an assumption, new ideas to 
explain the discrepancy could well be in order. Finally, we 
have demonstrated that our conclusions are robust when we 
exclude data with Galactic latitudes |6| < 30, which indi- 
cates that the information (or lack of information) in the 



Galactic plane is not responsible for the missing large scale 
power. 
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